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Abstract 

In this letter we study new relevant phenomenological consequences of the right-handed heavy neutrinos 
with masses at the 0(1) TeV energy scale, working within the context of the Inverse Seesaw Model that 
includes three pairs of quasi-degenerate pseudo-Dirac heavy neutrinos. We propose a new exotic signal 
of these heavy neutrinos at the CERN Large Hadron Collider containing a muon, a tau lepton, and 
two jets in the final state, which is based on the interesting fact that this model can incorporate large 
Lepton Flavor Violation for specific choices of the relevant parameters, particularly, the neutrino Yukawa 
couplings. We will show here that an observable number of prjj exotic events, without missing energy, 
can be produced at this ongoing run of the LHC. 


1 Introduction 

One of the most interesting properties of the Inverse 
Seesaw Model (ISS) is that it allows for the existence 
of heavy Majorana neutrinos with masses within the 
reach of the CERN Large Hadron Collider (LHC), 
and with couplings to the Standard Model (SM) par¬ 
ticles that can be sizable, therefore leading to ob¬ 
servable production and decay rates. The version of 
the ISS that incorporates three pairs of SM-singlet 
heavy neutrinos, one pair per generation, is espe¬ 
cially appealing because, on the one hand, it can 
explain the observed small neutrino masses and neu¬ 
trino oscillation data and, on the other hand in con¬ 
nection with the LHC physics, it can lead to new 
interesting phenomenology with experimental signa¬ 
tures distinct from those of the SM particles. The 
most frequently studied signatures of heavy neutri¬ 
nos are those related to their Majorana nature [1-3] 
and, in particular, the most characteristic signal is 
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the same-sign dilepton plus two jets events which is 
being searched for at the LHC. The rates of this type 
of events within the ISS are directly related to the 
Majorana mass matrix px that introduces the Lep¬ 
ton Number Violation in the model, which in turn is 
built in as to precisely fit the observed light neutrino 
masses via a low scale realization of the Seesaw Mech¬ 
anism where Mn g ht — ninM^ 1 pxM^m^ and mu 
and Mr are the Dirac and extended right handed 
(RH) neutrino sector mass matrices respectively. The 
usually assumed hierarchy among the three input 
scales, px <C wd "C -Mr, then produces the inter¬ 
esting mentioned features of the ISS. Furthermore, 
due to the smallness of px , the heavy neutrinos com¬ 
bine in nearly degenerate pairs to form pseudo-Dirac 
fermions with masses close to Mr and with a split¬ 
ting of order O(px)- This pseudo-Dirac character of 
the heavy neutrinos has also been explored in the lit¬ 
erature in connection with the appearance of other 
interesting signals at the LHC, like the trilepton final 
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state [4,5] , and other multi-lepton signals [ 6 ] . 

We propose here a new exotic signal of the ISS 
heavy neutrinos at the LHC that is based on another 
interesting feature of the ISS, the fact that it can in¬ 
corporate large Lepton Flavor Violation (LFV) for 
specific choices of the model parameters, particularly 
the Yukawa couplings Y v , and lead to a new rich phe¬ 
nomenology, including the appearance of rare pro¬ 
cesses like the radiative decays, k —> Ij'y (i ^ j) and 
others. Moreover, it has been shown recently [7] that 
sizable branching ratios for the LFV Higgs boson de¬ 
cays, like H —» fir, can also be produced within this 
ISS context, that can be enhanced in its supersym¬ 
metric version [ 8 ] for specific choices of the model pa¬ 
rameters, reaching values which are compatible with 
the present excess observed by CMS and ATLAS in 
H —>• [it [9,10], while respecting the present upper 
bounds on the related radiative decay r —> /ry [ 11 ]. 
Our specific proposal here is to look at rare LHC 
events with one muon, one tau lepton and two jets in 
the final state that are produced in these ISS scenar¬ 
ios with large LFV, and that presumably will have 
a very small SM background. This letter then sum¬ 
marizes our computation of the rates for these exotic 
firjj events due to the production and decays of the 
heavy quasi-Dirac neutrinos at the LHC within the 
ISS. 


2 ISS scenarios with large LFV 


We consider the ISS model that extends the SM with 
two kinds of fermionic gauge singlets: RH neutrinos, 
denoted here by Vr and coupled to the vr via Yukawa 
interactions, and X, with opposite lepton number. 
We assume one pair of them per generation. The 
relevant ISS Lagrangian for the mass terms is: 


Ass = A'*',,, -MpgXj - ',/VYf.Y, • h.c., 

( 1 ) 

where L is the SM lepton doublet, $ is the SM Higgs 
doublet with (4>) = v = 174 GeV, $ = ia 2 &*, Y v is 
the 3x3 neutrino Yukawa coupling matrix, Mr is a 
lepton number conserving complex 3x3 mass matrix, 
and jJx is a Majorana complex 3x3 symmetric mass 
matrix that violates lepton number conservation by 
two units. The diagonalization of the relevant 9x9 
neutrino mass matrix, containing the Dirac mass ma¬ 
trix, mo = Y u (<&), and Mr, 
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provides the 9 physical neutrino states rij(j = 1, ..9): 
three light neutrinos vi, 2,3 that are identified with the 


ones experimentally observed, and six heavy neutri¬ 
nos W ,..6 which are the object of our interest here. In 
order to match the predictions of this ISS model with 
the low energy neutrino data, one considers small fix 
given by the following parameterization in terms of 
the light neutrino mixing matrix C/pmns [12]: 

l l x = Mr m R Dp MNS m,,f7p MNS ml Mr. (3) 


This parametrization, which works quite well in the 
fix rriR <C Mr limit, was introduced in ref. [7] 
and allows us to use Y v and Mr as input parameters 
of the model. For the rest of this letter, we focus 
on the simplest case with diagonal Mr matrix and 
degenerate entries, Mr, 2 3 = Mr. Within this men¬ 
tioned small fix limit, the heavy neutrinos can be 
grouped into pseudo-Dirac pairs with nearly degen¬ 
erate masses. Both the light neutrino masses and the 
small differences among the quasi-degenerate heavy 
neutrino masses are governed by the small size of 
fix, basically following the same pattern as in the 
one generation case where: 
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For the present work, the corresponding pairings of 
the six heavy neutrinos are given by: N- t i-i, IV3/4 and 
A 5 / 6 , with m Nl/2 < TOat 3/4 < m N5/6 . Their corre¬ 
sponding masses in the Mr limit are close to 

Mr which will be taken here in the range accessible 
to the LHC, namely < 0(1 TeV). 

Regarding the heavy neutrino interactions, we 
give here just the relevant ones, namely, those to W 
gauge bosons, which enter both in the heavy neu¬ 
trino production and in the decays of our interest, 
and those to Z gauge and Higgs bosons, which enter 
in the computation of the total heavy neutrino width: 
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where, assuming the charged leptons to be diagonal, 
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and U 1 ' is the 9x9 neutrino rotation matrix 
that diagonalizes the mass matrix ( 2 ) according to 
U L ' T MissU v = diag(m ni ,... ,m„ 9 ). 
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Since we are interested in the ISS with large LFV, 
and more specifically producing large r — /x transi¬ 
tions, we will work with some benchmark scenarios 
that produce these wanted features. In particular 
we will use the three benchmark scenarios proposed 
in [7] that are defined by the following input Yukawa 
matrices, Y v \ 
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These scenarios have been explored fully in [7] 
and display large LFV ratios for processes involv¬ 
ing r — /x transitions (this motivates the given name, 
Y v = Y^), whereas they strongly suppress pro¬ 
cesses involving r — e and /x — e transitions, which can 
be easily understood with the use of the approximate 
formulas presented also in [7]: 
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(13) 

For / values close to its perturbativity limit, 

/ = \/67r, and M R ~ 0(1 TeV) these lead 

to rates of up to BR(iJ —► jif) ~ 10~ 5 while 
respecting the present bounds on the radiative 
decays, BR(/x —► ey) < 5.7 x 10"~ 13 [13], 
BR(r —► /xy) < 4.4 x 10 -8 and BR(r —>■ ey) < 

3.3 x 10 -8 [11]. Moreover, in the case of the SUSY 
version of this ISS model, the Higgs decay rates in¬ 
duced by the SUSY loop contributions have been 
found in [8] to be as large as the present LHC sensi¬ 
tivity of about BR (H —> /xr) ~ 10 -2 [9,10]. 

As we see from eq. (12), the constrains from LFV 
radiative decays set upper bounds on the combina¬ 
tion \(Y v Yj)ij\, with i ^ j. In terms of the heavy 
neutrino flavor mixing matrix, H/at, they constrain 
the combination \{B lN B' lN )ij\, but not the B l N . it¬ 
self, which is the relevant element that controls the 
flavor pattern of each of the heavy neutrinos. 


Figure 1: Heavy neutrino flavor mixing, according to 
eq.(14), for the three benchmark scenarios in (11). 


In figure 1 we show the flavor content of each 
heavy neutrino for the three benchmark scenarios in 
(11), where the length of the colored bars is calcu¬ 
lated as 


SlNi = 


\BlNi \ 2 

i B iNi \ 2 


(14) 


and, therefore, represents the relative mixing of the 
heavy neutrino V,; with a given flavor l. We learn 
from fig.l that, although the three scenarios share the 
property of suppressing the LFV /x-e and r-e rates 
while maximizing the r-/x ones, the heavy neutrino 
flavor mixing pattern is different in each scenario. 
We also see that some heavy neutrinos carry both 
/x and r flavors, specially in the first two scenarios, 
pointing towards signals with both /x and r leptons. 
It should be further noticed, that the values of these 
Bin mixing parameters are determined within the 
ISS in terms of the input mr> and Mr mass matrices 
and, therefore, in the limit me < Mr they are sup¬ 
pressed as Bin ~ 0(to£>M^ ). In addition, they are 
constrained from ElectroWeak Precision Observables 
(EWPO) [14] by: 


\B eN \ 2 < 3.0 x 10“ 3 , 

(15) 

|^at| 2 < 3.2 x 10- 3 , 

(16) 

| B r n | 2 < 6.2 x 10" 3 . 

(17) 
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3 Numerical Results 


Heavy neutrinos with masses of the TeV order can 
be produced at present and future colliders, in par¬ 
ticular in the second run of the LHC. The dominant 
production mechanism in this case is the Drell-Yan 
(DY) process, fig-(2) left, where the heavy neutrino is 
produced in association with a charged lepton. The 
7 W fusion, fig. ( 2 ) right, also produces the same sig¬ 
nal with two extra jets and, in fact, can be relevant 
especially for large neutrino masses in the 0(1 TeV) 
energy range [15,16]. 

In order to estimate the heavy neutrino produc¬ 
tion at the LHC, we have implemented the model 
in MadGraph 5 [17] using FeynRules [18]. Follow¬ 
ing ref. [16], we have used a AT-factor of 1.2 for 
the DY-process and split the 'yW process into three 
regimes characterized by the virtuality of the pho¬ 
ton 1 : elastic, inelastic and deep inelastic scatter¬ 
ing (DIS) regimes. In particular, we have set the 
boundaries between these three regimes to A® las = 
1.22 GeV and Aj^ IS = 15 GeV. In order to regularize 
possible collinear singularities, we have also imposed 
the following cuts to the transverse momentum and 
pseudorapidity of the outgoing leptons: 


p l T > 10 GeV, I??*] <2.4. (18) 


The results for the three benchmark scenarios in 
(11) are shown in fig. 3, where the dominant DY pro¬ 
duction cross sections normalized by / 2 are plotted 
as functions of the mass parameter Mr. We see that 
the production cross sections can be of the fb or¬ 
der, reachable at the LHC, for masses Mr < 600 
GeV. Notice that the results are always equal for the 
pseudo-Dirac pairs, since their Majorana character 
plays a subleading role in their production. 



Figure 2: The two main processes, Drell-Yan and r )W 
fusion, producing exotic events via heavy neu¬ 
trino production and decay at the LHC. 
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Figure 3: Heavy neutrino DY-production at the LHC 
for the three scenarios. The cross sections are nor¬ 
malized by the Yukawa coupling intensity / in (11). 
Processes not shown are negligible. 

We can also learn that the flavor of the associated 
charged lepton is different depending on the heavy 
neutrino produced and the scenario considered, and 
that this pattern can be understood looking at the 
mixing in fig.l. Thus, IV 3/4 are mainly electronic 
neutrinos in the first and third scenarios, therefore, 
they are basically produced exclusively with elec¬ 
trons. N 1/2 are equally produced with muons and 
taus in the first scenario, dominantly produced with 
electrons in the second scenario, and mainly produced 


1 We warmly thank Richard Ruiz, Tao Han and Daniel Alva for their generous help and clarifications in the implementation 

of the 7 W process. 


4 


























with taus in the third scenario. N 5 / e are equally pro¬ 
duced with muons and taus in the first and second 
scenarios, and mainly produced with muons in the 
third scenario. 

Once the heavy neutrinos are produced, they will 
decay inside the detector. As mentioned above, in 
the limit Mr > mjj the heavy neutrino masses are 
close to Mr, with small differences of 0{m 2 D M^ j 1 ) be¬ 
tween the different pseudo-Dirac pairs and, therefore, 
assuming they are practically degenerate, their decay 
into each other should be suppressed, with the dom¬ 
inant channels being Nj —► Zvi,Hvi,W ± lf. The 
relevant decays into W + l~ or W~lf have a partial 
width given by: 



M r (GeV) 



We have also seen that, in agreement with [3,4], in the 
large Mj limit, the heavy neutrinos decay with equal 
branching ratios to all the bosons in all the scenarios 
if we sum over all the flavors, namely, 

BR (Nj -> hv) = BR (Nj -> Zv) = 

BR (Nj -> W+l~) = BR {Nj ->W~l + ) = 25% . 

( 21 ) 

It is, however, interesting to study the rich flavor 
structure of the decay products, which depends on 
the decaying heavy neutrino and the scenario we are 
considering. Like in the production, the flavor pref¬ 
erence of the decays to W ± P, which are the ones 
relevant to this study and are given in fig. 2 , also fol¬ 
lows the same pattern as in fig. 1 . Therefore, we can 
expect the production and decay of the heavy neu¬ 
trinos to lead to exotic /irjj events with no missing 
energy and Mjj ~ M\y, with Mjj the invariant mass 
of the two jets. 

Using the narrow width approximation, the total 
cross section of the exotic events we are interested in 
is given by: 

cr(pp ->• nrjj) = 

6 

^ jer(pp —y N i p ± )BR(N i —► W ± t t ) 

i=1 


/W R (GeV) 



M r (GeV) 


Figure 4: Number of exotic /irjj events at the LHC 
for the three scenarios and for three values of /. 
For each /, the bottom solid line is the prediction 
of /irjj events from DY and the upper lines on top 
of each of the three shadowed regions are the pre¬ 
dictions after adding the /irjj events from 7 W, for 
p™ ax = 10 , 20,40 GeV from bottom to top, respec¬ 
tively. The p™ ax cut is only applied to the two extra 
jets in the events from 7 W. The upper red shadowed 
areas are excluded by r —> /ry. 


+ a(pp -> N i T ± )BR(N i -> W r± /t =F )} 
xBR (W ± ^jj). (22) 


Figure 4 shows the expected number of exotic 
events /irjj at the LHC for an integrated luminosity 
of C = 300 fb ” 1 at -y/s = 14 TeV. The lower solid 
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lines for each choice of / are the number of events 
considering only the DY-production. 

Moreover, 7 W fusion processes can also con¬ 
tribute to this kind of exotic events if the pr of the ex¬ 
tra jets are below a maximum value p™ ax and, there¬ 
fore, they can be considered as soft or collinear jets 
which can escape detection. In this case the predicted 
total number of exotic events are the sum of the 
events produced by DY and 7 JU channels. These to¬ 
tal contributions for different values of p^ ax = 10 , 20 
and 40 GeV are shown as the border lines on top of 
the shadowed areas with gradual decreasing intensity 
above each solid line. In addition we have included 
in the plots red shadowed areas that represent the re¬ 
gions excluded by the experimental upper bound on 
BR(r —> fi'y). We can see that, after considering all 
the constraints, the three scenarios lead to an inter¬ 
esting number of 0 ( 10 — 100 ) total prjj exotic events 
for the range of Mr studied here of [200 GeV, 1 TeV]. 

The SM backgrounds for events with two leptons 
of different flavor have been studied in ref. [19]. How¬ 
ever, a high efficiency in the r-tagging and a good 
reconstruction of the W boson invariant mass from 
the two leading jets would help in reducing the back¬ 
ground. In that case, the main background would 
come from processes with photons or jets misidenti- 
fied as leptons, mainly from W/Z + 7 *, W/Z + jets 
and multijet events with at least four jets with one 
of them misidentified as a muon and another one as 
a tau; and from Z/"f* + jets —> p + p~ + jets if one of 
the muons is misidentified as a r candidate. Never¬ 
theless, a dedicated background study for these par¬ 
ticular prjj exotic events is beyond the scope of this 
letter and it will be done in a future work. 

4 Conclusions 

In this letter we have proposed a new interesting way 
to study the production and decay of the heavy neu¬ 
trinos of the ISS in connection with LFV. We have 
presented the computation of the predicted number 
of exotic jiTjj events which can be produced in these 
ISS scenarios with large LFV by the production of 
heavy pseudo-Dirac neutrinos together with a lepton 
of flavor £, both via DY and 7 W fusion processes, 
and their subsequent decay into W plus a lepton 
of different flavor. We have concluded that, for the 
three benchmark scenarios studied here, a number of 
0 ( 10 — 100 ) total prjj exotic events without missing 
energy can be produced at this ongoing run of the 
LHC, for values of Mr from 200 GeV to 1 TeV. Sim¬ 
ilarly, other rare processes like rejj or pejj could 
be produced within other ISS scenarios with large 
LFV, although for the latter ones, while being inter¬ 
esting since they could provide in addition observable 


CP asymmetries [20] , the number of events would be 
strongly limited by the p —► ey upper bound. Of 
course, a more realistic study of these exotic events, 
including detector simulation, together with a full 
background study should be done in order to reach a 
definitive conclusion, but this is beyond the scope of 
this letter and will be addressed in a future work. 
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